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Abstract 

In this paper we argue that having available higher order fully exclusive Monte Carlo 
generator for Higgs boson production from heavy quark fusion will be mandatory for data 
analysis at LHC. The H — > rr channel, a key for early discovery of the Higgs boson in 
the MSSM scenario, is discussed. With simplified example and for ran = 120 GeV we 
show, that depending on choice among presently available approaches, used for simulation 
of Higgs boson production from bbH Yukawa coupling, final acceptance for the signal 
events being reconstructed inside mass window may differ by a factor of 3. The spread is 
even larger (up to a factor of 10) for other production mechanisms (promising for some 
regions of the MSSM parameter space). The complete analysis, which necessarily will 
add stringent requirements for background rejection (such as identification of b-jet or veto 
on b-jet) and which will require statistical combination of samples selected with different 
selection criteria may only enhance the uncertainty. 
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1 Introduction 



The search for the Higgs boson is one of the primary task of the experiments at LHC 
[1,2]. Several studies discussed there, concluded both for the Standard Model and Minimal 
Supersymmetric Standard Model, that Higgs boson can be discovered at LHC, with the 
high margin of significance, for the full interesting range of its mass. 

One of the key signature for Higgs boson discovery, is with its decay into pair of r- 
leptons. For the MSSM Higgs scenario, decay into r-pair is strongly enhanced at large 
tan/3, the parameter of the supersymmetric models. The promising (it depends on the 
Higgs boson mass and tan (3) production mechanisms are: the Higgs boson production in 
association with the b-quarks, inclusive gluon-gluon fusion and, for the lightest MSSM 
Higgs, also vector-boson-fusion. 

The analyses prepared for the MSSM Higgs boson searches in ATLAS, presented in 
[1], rely mostly on the lepton-hadron decays of the r-pair, that is the case, when one of 
the r-leptons, from the Higgs boson decay, subsequently decays leptonically and another 
one hadronically (hadron-hadron case was discussed only for the Higgs boson of large 
mass). The lepton-lepton channel contribution, to the Higgs boson discovery potential, 
was discussed in [2], in context of the CMS detector. 

On theoretical side, the inclusive cross-section at large tan/3 is dominated by the 
bottom-quark fusion process bb — > H. For the first time it was discussed in [3]. Recently, 
important progress has been achieved. The total cross-section for the bb — > H process has 
been evaluated to the next-to-next-to-leading order (NNLO) [4] in the so called variable 
flavour number scheme (VFS) [3, 5, 6]. The result of this NNLO calculation shows almost 
no scale dependence. The inclusive bb — > H cross-section was obtained also at the next- 
to- leading (NLO), fixed order calculations for the parton level process gg,qq — ► bbH [7, 8]. 
There, the fixed flavour number scheme (FFS) was used. 

Results obtained in these two schemes seem to be compatible now with each other, and 
show that there is actually no large difference between the NLO fixed order results and use 
of the b-quark structure functions, if the proper factorisation scale for this process is used 
[6]. Without the NNLO calculations [4], understanding consistencies between FFS and 
VFS approaches would be on much weaker foundation. This statement concludes what 
was discussed since a long time [5, 6, 9, 10]. The fixed order calculation shows a substantial 
scale dependence and sufficient control of the residual large uncertainties, available only 
at higher orders, was necessary to compare results of the two approaches. 

One should be well aware, that for a given decay channel of r-pair, the statistical 
combination of the discovery evidence for events with one identified 6-jet and events with 
no identified 6-jet at all will be required, to achieve maximal (sufficient) sensitivity to 
the Higgs boson. This defines, what are the relevant hard processes for the analyses, as 
designed in [1]. For example if the identified final state includes single bottom quark (jet), 
then the relevant lowest order hard process in the VFS approach should be gb — > bH [5]. 
The cross-section for the gb — ► bH production has been also computed at NLO [9] and the 
residual uncertainties due to the higher order corrections are small. In the FFS scheme 
the relevant lowest order hard process will still be the gg,qq — > bbH. 

At early stage of LHC operation, also statistical combination of signatures with dif- 
ferent r decay channels will need to be performed. In addition, in the MSSM scenario, 
signals from different Higgs bosons (h, H and A) not degenerated in mass, might nonethe- 
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less overlap in the same mass window. 

Although the NLO and even NNLO calculations became available for the integrated 
cross-sections, and impressive progress has been achieved in understanding consistencies 
between (VFS) and (FFS) approaches [4, 5, 6, 9, 10], and even though some results on 
transverse Higgs momentum distribution are also discussed in the literature (see. e.g. [11]), 
only the LO matrix element + parton shower approach is available for the full event 
generation which allow for the subsequent detector simulation 1 . 

As we will attempt to illustrate in this paper, the foreseen by experiments signal recon- 
struction procedures are very sensitive to the topology of the signal production process. 
Therefore, good understanding, from the theoretical perspective, of the Higgs production 
topologies will be mandatory. This can be achieved only with fully exclusive higher order 
Monte Carlo generator, what we hope will be visible from the present paper. 

In this paper, as an example we will use the SM-like 120 GeV mass Higgs boson. We 
will compare reconstruction efficiencies and final resolutions for the different hard processes 
and for lepton-lepton and lepton-hadron r decay modes. 

Our paper is organized as follows. In section 2 we discuss Higgs boson production using 
different hard processes and the appropriate cross sections. In sections 3 and 4 we review 
basic reconstruction and selection properties to be used by experiments. This in particular 
will explain how idealized signatures translate into more realistic ones. In section 5 we 
collect numerical results for the {I I pj: lss ) signature, originating from the case when both 
r-leptons decay leptonically. Similarly, in section 6 we collect numerical results for the 
(£ r-jet p™ lss ) signature, originating from the case when one r decays leptonically and 
another one hadronically. Finally, conclusions, section 7 close the paper. 

2 Different hard processes 

Let us start discussion by presenting the naive Table 1 with cross-sections calculated for 
three different hard processes: the (2 — ► 1) process bb — > H, the (2 — ► 2) process gb — ► Hb 
and the (2 — ► 3) process gg, qq — > bbH, as obtained from PYTHIA 6 . 2 [12] generator, 
according to its default initialisations 2 . 

Factor 4 difference in normalisation can be reported between bb — > H and gg, qq — > bbH 
hard processes. The first one represents lowest order term in (VFS) scheme, second the 
lowest order term in the (FFS) scheme for the Higgs boson production mediated by the 
bbH Yukawa coupling, which does not rely on the 6-quark PDF's. Given recent clarification 
in [6] this difference could be minimised by using the proper factorisation scale. 

The question, whether differential distributions can impose significant effects on how 
in experimental conditions, the Higgs boson signature is expected to be defined, was not 
adressed so far in theoretical calculations. The experimental studies have shown clearly 

1 Here we refer explicitly to PYTHIA or HERWIG generators, which provide generation of bb — > H, gb — ► bH 
and gg — ► bbH in the lowest order of the hard process only. We are not aware of any implementation in 
form of the Monte Carlo event generator of the complete calculations of [4]. 

2 These means choices for the QCD factorisation scale, minimum bias model, parameters of the shower 
evolution, etc. etc. We think that, for example, discussion of consequences of the possible alternative 
choices is beyond the scope of this paper, it could only dilute the aim of the paper. The appropriate 
theoretical framework, necessary for such a discussion need to be established first. For the structure 
functions we have used CTEQ5L parametrisation. 
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that the impact of the topological features of the production process can be significant on 
the overall efficiency of standard reconstruction procedure, see ref. [1] pages: 746, 747. 

One of the characteristic for event topology differential distribution is the transverse 
momenta of the Higgs boson, as it determines directly the average transverse momenta 
and angular separation of the decay products. In Fig. 1 we show Higgs boson transverse 
momenta distribution p^" s as generated with different hard processes used. The average 
< p^ l " s > distribution for bb — > H process is 23 GeV, for the gb — ► bH process is 31 GeV, 
and finally for the gg — ► bbH process is 28 GeV. The differences in event topologies seems, 
at a face value, not to be dramatic but we will study nonetheless their impact further in 
the paper. 



Hard Process 


ax BR [fb] 


bb H(-> tt) 


7.2 


gb bH(-+ tt) 


4.1 


99, qq bbH{-^ tt) 


1.7 



Table 1: Cross-section for signal production with b-quark Yukawa coupling. Results for 
there different hard processes are collected. Branching ratios of H — > tt is included. The 
cross section is for 120 GeV mass SM Higgs boson at 14 TeV pp collision, as generated 
with default initialisations in PYTHIA. 




Figure 1: Thep T % " s distribution for the three hard processes used: bb — > H (left), gb — ► bH 
(middle) and gg — ► bbH (right). Distribution normalised to total a x BR (pb). Each of 
these distribution is expected to include significant deficiencies. In the first case generation 
of pt above certain value (max. shower scale) is not performed, the later two, are valid 
only when respectively one/two final state b-quarks havepT above Sudakov turnover. These 
condition do not translate into constraint on p^ l " s in a straightforwad manner. 

To perform realistic studies it was not enough to generate Higgs production process. 
For that purpose, PYTHIA [12] Monte Carlo was used. Essential element consisted 
of simulating the decays of the r-leptons. For that purpose TAUOLA [13, 14, 15] and 
PHOTOS [16, 17] Monte Carlos were used. The three segments of the generation were com- 
bined using interface presented in [18]. In this way, correct simulation of the r-decays and 
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the radiative photons emission from the Higgs boson decay cascade, relevant in particular 
for the hadron-lepton mode were achieved. Also full spin correlation effects in the Higgs 
boson decays were assured. 

3 Reconstruction of the basic experimental signatures 

Reconstruction of the basic experimental signatures: leptons, r-jets, missing transverse en- 
ergy, have been done with the fast simulation of the simplified LHC detector, AcerDET [19]. 
Although to large extend it is representing the best performance of the detector, we believe 
that it is fairly adequate for the studies presented in this paper. The principle of AcerDET 
operation is quite similar to the official fast simulation package of ATLAS Collaboration, 
ATLFAST [20]. In contrary to the latter it relies on only those parameters of the detector 
which are available in the literature, eg [1], thus public. It is however missing regular 
internal updates of ATLFAST. 

The key experimental issues, which are relevant for discussed analyses can be shortly 
summarised as: • Both discussed channels have leptons (e or (i) in the final state and 
can be triggered by either single- or the dilepton trigger. • The hadronically decaying 
tau's can be identified in the detector with good efficiency and purity. The important 
ingredient in the r-jet identification is the profile of the energy deposition in calorimeter 
and the number of tracks pointing to the calorimeter cluster. • Missing transverse energy, 
p™ ss , is calculated after summing up all reconstructed visible signatures in the final states, 
including also cells and clusters in the calorimeter. 

We have checked that if the ATLFAST is used in simulation instead of AcerDET, a bit 
worse mass resolution and efficiency for the r-jet reconstruction is obtained. This does 
not lead to changes in final conclusions of the paper. For more detailed discussion on 
the comparison of full ATLAS detector simulation and AcerDET results, relevant for the 
present analysis, see [19]. 

4 Selection criteria 

Discussed here selection criteria are the minimal one for triggering of a given event (at 
least one lepton in the final state) and for a good reconstruction of the r-pair mass. 
Reconstruction of the invariant mass of the r-lepton system is made under assumption 
that the r-lepton is massless and decays collinearly 3 [21, 1, 2]. The procedure which have 
been used in [1, 2] is equivalent, for the physical solutions, to the one used in [22, 23] 
and described below. Fractions of the two t's momenta which are carried by measured 
visible decay products, x T1 , x T2 , can be calculated from solving equations of conservation 
of the transverse momenta components of the Higgs boson, reconstructed from the visible 
products (leptons, r-jets) and missing transverse energy. The physical solutions (events 
with resolved neutrinos) are those for which < x T1 ,x T . 2 < 1. For events with physical 
solution, the invariant mass of the system of the visible decay products of r's is calculated 
and the invariant mass of the r-system is expressed as m TT = m v i s /^/x T1 ■ x T2 . 

3 We have checked that this assumption leads to about 2.5-3.5 GeV contribution to the gaussian width 
of reconstructed invariant mass of the rr system. 
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We start with selection procedure as denned in [1] and refined in [24] for the lepton- 
lepton final state. The basic selection consists of requiring: two isolated leptons within 
acceptance of the detector; threshold on the minimal angular separation between leptons; 
threshold on the reconstructed missing transverse energy and kinematics of events which 
allow to resolve equations for neutrinos momenta, i.e. for invariant mass of the r-pair 
system to be calculable. The additional selection, also introduced already in [1, 24], is 
necessary to optimise resolution of the reconstructed invariant mass of the r-lepton pair. 

Slightly modified selection is used for lepton-hadron final state, the basic selection 
criteria have to be adjusted, already from the beginning, to reject expected reducible 
backgrounds (ti, W+j). 

For the explicit definition of basic selection and additional selection cuts 4 see later: 
Tables 2 and 5 respectively in sections 5 and 6. 

5 The (£ I P rp tss ) final state 

In Fig. 2 we show distributions of the kinematical variables in the transverse plane: 
sin(A^), p™ ss , cos(A0«), used for event selection. One can observe some differ- 
ences in the shape of the distributions, which do not seem very significant (except for 
cos(A^)), but nevertheless lead to the cumulated effect on the acceptances of order of a 
factor 3 (see Table 2) for signal events in lepton-lepton mode, generated with different LO 
processes. The very initial acceptances are comparable (first line of Table 2), divergences 
start with adding selection and reconstruction requirements. 



Selection 


bb^H 


gb^bH 


bbH 


99 


qqH 


Basic selection 












2 iso £, p l T > 15 GeV 


18.6 • 10~ 2 


18.4- 10~ 2 


19.2 • 10~ 2 


18.6 • 10~ 2 


21.3 • 10~ 2 


|sm(A<^)| > 0.2 


9.3 • 10~ 2 


10.1 • 10~ 2 


10.1 • 10~ 2 


10.4 • 10~ 2 


19.0 • 10~ 2 


p^ lss > 15 GeV 


5.5 • 10" 2 


7.0 • 10" 2 


6.4 • 10" 2 


7.9 • 10" 2 


17.7- 10~ 2 


resolved neutrinos 


4.4 • 10~ 2 


5.8 • 10~ 2 


5.1 • 10~ 2 


7.0 • 10~ 2 


16.5 • 10" 2 


Additional selection 












pip tss > 30 GeV 


1.4- 10~ 2 


3.0 • 10~ 2 


2.2 • 10~ 2 


4.2 • 10~ 2 


13.4 • 10~ 2 


cos(A(p u ) > -0.9 


1.0 • 10~ 2 


2.5 • 10" 2 


1.8 • 10~ 2 


3.8 • 10" 2 


12.8 • 10~ 2 


Ru < 2.8 


9.5 • 10~ a 


2.5 • 10~ 2 


1.8 • 10~ 2 


3.7- 10" 2 


12.7- 10~ 2 



Table 2: The cumulative acceptances of the selection criteria for different approaches of 
modeling production process. For each subsequent line effect of the additional cut off is 
added. 

From the Table 2 we can read, that a large systematic theoretical uncertainty must be 
assumed for the efficiency of the selection and reconstruction procedure, as a consequence 

4 Angle <f> is measured in the plane transverse to the beam axis. The Ru, Re T -jet denote cone separation, 
expressed by the difference in pseudorapidity f] and angle <j>, R = y (Acj>) 2 + (Arj) 2 . The transverse mass 

in lepton-hadron channel is calculated as mt, m isa = \J (Pt) 2 + (p™ ISS ) 2 - For more comments on definition 
of variables see [24] . 
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of the different choices in modeling of the production process topology. For illustration, 
we include results for the other production processes 5 : gluon- fusion gg — ► H, and vector- 
boson-fusion, qq — > qqH as well. These processes will be included in the complete analysis 
of the Higgs signal observability as the channels are expected to be promising for some 
regions of the MSSM parameter space [22, 23]. The final spread on the efficiency is then 
even larger. 




20 40 60 80 100 120 140 160 180 200 

pf s (GeV) 

]4 Mean =-0.6721 

' RMS =0.3087 



100 120 140 160 180 200 

pf(GeV) 

Mean =4.2301 
RMS =0.6108 



lax, nHlr^ 



20 40 68 80 100 120 140 160 180 200 

pf s (GeV) 

Mean = -0.404I 
RMS =0.5193 



,rK,. 



C0S(Af)(l,l)) 



COS(Ao(U! 



COS(A0)(l,i; 



Figure 2: The characteristic kinematical distributions before respective selection, for differ- 
ent hard processes used: bb — > H (left column), gb — ► bH (middle column) and qq — ► bbH 
(right column). The shaded (yellow) area will be accepted by selection. Distributions nor- 
malised to total a x BR [pbj. Every line of the plots correspond to variables as used in 
cuts: lines 2, 3 and 6 of table 2. Symmetric shape of the distributions of the first line of 
the plots provides technical test of the simulation. 

5 For the production processes gg — » H and qqH, which we will also include later in our discussion of 
acceptance, default initialisations of PYTHIA 6 . 2 are used as well. 
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Selection 


bb^H 


gb^bH 


bbH 


99 


qqH 


Basic selection 


17.3 GeV 


14.1 GeV 


16.0 GeV 


11.6 GeV 


8.8 GeV 




(63.6 %) 


(68.6 %) 


(67.4%) 


(79.1%) 


(88.7%) 


p^ lss > 30 GeV 


15.8 GeV 


11.9 GeV 


13.0 GeV 


9.8 GeV 


8.3 GeV 




(73.0 %) 


( 77.7 %) 


(76.9%) 


(86.5%) 


(91.3%) 


cos{A(j) U ) > -0.9 


14.5 GeV 


11.4 GeV 


12.3 GeV 


9.3 GeV 


8.1 GeV 




(80.8 %) 


( 82.6 %) 


(82.0%) 


(90.0%) 


(93.1%) 


Ru < 2.8 


14.4 GeV 


11.4 GeV 


12.1 GeV 


9.2 GeV 


8.1 GeV 




(81.6 %) 


( 83.2 %) 


(82.9%) 


(90.6%) 


(93.4%) 



Table 3: Resolution of the Gaussian fit to the reconstructed invariant mass of the tt system 
for different approaches of modeling production process. In brackets is shown acceptance 
within mass window of mjj ± 20 Ge V. 




Figure 3: The reconstructed mass of the t's pair, for the gg — > bbH events. Distribution 
after cos(A</>«) selection cut off is shown. Only those events which are reconstructed in 
constrained mass window, eg. mu ± 20 GeV, will be considered as contributing to the 
Higgs signature. 

In Table 3 we compare resolution of the Gaussian fit to the reconstructed invariant mass 
of the r-pair system, and acceptances inside mass window of mu ± 20 GeV. The natural 
width is negligible for the Higgs boson with 120 GeV mass, the non-zero resolution spread 
comes exclusively from the reconstruction procedure (assumptions of the collinearity in 
the r's decay, and p™ ss reconstruction). As an example, in Fig. 3 we show reconstructed 
invariant mass for the gg —>■ bbH events. Resolutions for bb —>■ H and gb — > bH topologies 
(Table 3) differ up to 30%. This additional effect will enhance the impact of choice of the 
hard process on the total acceptance, on top of the one, introduced by the differences in 
the selection efficiencies discussed in Table 2. 
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Estimates for the total cumulative acceptance inside mass window, as presented in 
Table 4, give clear indication of the size of the systematic uncertainty which should be (at 
present) assigned to the predictions on the expected number of signal events. 

Note that the very initial acceptances (first line(s) in Table 2) were comparable for 
all discussed processes. Then the selection was designed for signal reconstruction only. 
Additional features of topologies, like presence of extra hard partons was not explored yet. 
The difference for the total cumulated acceptances, which arised after additional selection 
necessary for the signal reconstruction inside mass window, is a factor of 3 for bb — > H 
versus gb — > bH topologies. It is even a factor of 4 for gg — > H versus bb — ► H topologies. 
The additional selection as it was already pointed out, is necessary for resolution of the 
reconstructed invariant mass of the r-lepton pair. 



Selection 


bb^H 


gb^bH 


bbH 


99 


qqH 


Basic selection 


2.8 • 10~ 2 


4.0 • 10~ 2 


3.4 • 10~ 2 


5.6 • 10- 2 


14.6 • 10~ 2 


p?p iss > 30 GeV 


1.0 • 10~ 2 


2.3 • 10~ 2 


1.7- 10~ 2 


3.7- 10~ 2 


12.2 • 10~ 2 


cos(A(f> u ) > -0.9 


8.2 • 10~ 3 


2.1 • 10~ 2 


1.5- 10~ 2 


3.4- 10~ 2 


12.0 • 10~ 2 


R u < 2.8 


7.8 • 10~ 3 


2.1 • 10~ 2 


1.5- 10~ 2 


3.4- 10~ 2 


11.9 • 10~ 2 



Table 4: The cumulative acceptance in the mass window mn ± 20 GeV, for different 
approaches of modeling production process. 

6 The I r-jet pfp tss final state 

Let us now turn to the second case, where one of the r's decays hadronically. We will 
proceed similarly as in the previous section. We start with Fig. 4 where we show distribu- 
tions of the kinematical variables used later for events selection: sin(A0^ T -jet), m i T miss ^ 
cos(A<^ r-jet)- One can observe some differences in the shape of these distributions, de- 
pending on the choice of the hard process (similarly as for lepton-lepton channel, the 
cos(A<^ T -jet) is the most sensitive to this choice). The cumulated effect on the accep- 
tances is of the order of a factor of 3 in this case (see Table 5) similarly as it was in the 
lepton-lepton channel. 

The Gaussian resolution for the reconstructed invariant mass of the rr system mass 
is specified in Table 6. Obtained resolution on average is 10% better than for the lepton- 
lepton channel, which might seem contrary to the intuition (we need to reconstruct r-jet, 
what is less precise than reconstruction of a lepton). The contradiction is due to the 
fact that the r-jet spectrum from r — > had v T decay is harder than the lepton spectrum 
from r — > tv T vt decay (in the second case we have two neutrinos). Thus, in the first case 
the assumption of the collinear decay for r lepton works better for the invariant mass 
reconstruction, simply because there is less of neutrino energy. 

Estimates for the total cumulative acceptance inside mass window, presented in Ta- 
ble 7, give clear indication of the size of the systematic uncertainty which should be as- 
signed to the predictions for the expected number of signal events. The difference for the 
total cumulated acceptance for the signal reconstruction inside mass window are similar 
in size to the one in the lepton-lepton final state. 
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sin(A4i(l,T-jet)) 



Mean 


= 15.67 


RMS 


= 15.17 



100 120 140 160 180 200 

m^ miss (GeV) 



Mean 


= -0.731: 


RMS 


= 0.2704 



sin(A$(l,T-jet)) 



Mean 


= 17.36 


RMS 


= 16.58 



100 120 140 160 180 200 

m^lGeV) 



Mean = -0.137' 
RMS =0.6574 



1 „ r^l 




0.16F 

0.14 

0.12; 



sin(At(l,t-jet)) 



Mean - 


16.39 


RMS = 


15.9 



20 40 60 80 100 120 140 160 180 200 

my iss (GeV) 



Mean =-0.335; 
RMS =0.5622 



cos(A«)(l,T-Jet)) 



cos(Ao)P,T-]el)) 



cos(A<t>(l,t-jet)) 



Figure 4: The characteristic kinematical distributions before respective selection, for dif- 
ferent hard process used: bb — > H (left column), gb — > bH (middle column) and qq — > bbH 
(right column). The shaded (yellow) area will be accepted by the selection. Distributions 
normalised to total a x BR [pb]. Every line of the plots correspond to variables as used in 
cuts: lines 2, 3 and 6 of table 5. Symmetric shape of the distributions of the first line of 
the plots provides technical test of the simulation. 

Let us stress that results presented in Sections 5 and 6 should be considered as il- 
lustration of the problem, and not as an optimised expected performance of the LHC 
experiments. 
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Selection 


bb^H 


gb^bH 


bbH 


99 


qqH 


Basic selection 












1 iso £, p e T > 20 GeV 
1 r-jet, p T T ~ jet > 30 GeV 


19.5 • 10~ 2 


19.3 • 10~ 2 


19.7- 10~ 2 


19.5 • 10" 2 


22.2 • 10~ 2 


\sin(A(f>e T -jet)\ > 0.2 


9.6 • 10~ 2 


10.5 • 10~ 2 


10.4 • 10~ 2 


10.9 • 10~ 2 


19.8 • 10~ 2 


rrij, < 50 GeV 


8.9 • 10" 2 


9.8- 10~ 2 


9.7- lO" 2 


10.1 • 10~ 2 


18.4 • 10~ 2 


resolved neutrinos 


5.5- 10~ 2 


6.6- 10~ 2 


6.0 • 10~ 2 


7.6 • 10" 2 


15.3- lO" 2 


Additional selection 










pip zss > 30 GeV 


9.1 • 10~ 3 


2.1 • 10~ 2 


1.4- 10 -a 


3.0- 10~ 2 


10.1 • 10~ 2 


cos(A4> e T _ jet ) > -0.9 


6.5 • 10" 3 


1.8 • 10~ 2 


1.2 • 10~ 2 


2.7- 10" 2 


9.8- 10" 2 


Re r-jet < 2.8 


6.1 • 10" 3 


1.7- 10~ 2 


1.1 • 10~ 2 


2.6- 10" 2 


9.7- 10~ 2 



Table 5: The cumulative acceptances of the selection criteria for different approaches of 
modeling production process. For each subsequent line effect of the additional cut off is 
added. 



Selection 


bb^ H 


gb -► bH 


bbH 


99 


qqH 


Basic selection 


18.3 GeV 


14.4 GeV 


16.8 GeV 


12.8 GeV 


10.0 GeV 




(62.2 %) 


(65.5 %) 


(66.0%) 


(71.5 %) 


(82.4 %) 


p^ lss > 30 GeV 


13.6 GeV 


11.4 GeV 


12.5 GeV 


9.8 GeV 


9.0 GeV 




( 62.5%) 


( 72.6 %) 


( 71.1%) 


( 78.2 %) 


(85.2 %) 


cos(A4> e T _ jet ) > -0.9 


13.1 GeV 


11.1 GeV 


11.9 GeV 


9.7 GeV 


8.9 GeV 


( 71.5%) 


( 79.0 %) 


( 77.4%) 


( 82.7 %) 


(87.1%) 


R(. r-jet < 2.8 


13.0 GeV 


11.1 GeV 


11.8 GeV 


9.7 GeV 


8.9 GeV 


( 73.4%) 


( 80.1 %) 


( 78.0%) 


(83.4 %) 


( 87.6 %) 



Table 6: Gaussian resolution for the reconstructed tt system for different approaches of 
modeling production process. In brackets acceptance within mass window ofmu ± 20 GeV 
is shown. 



Selection 


bb^H 


gb^bH 


bbH 


99 


qqH 


Basic selection 


3.4- 10~ 2 


4.3- 10~ 2 


4.0 • 10~ 2 


5.4- 10~ 2 


12.6 • 10~ 2 


tfP* 88 > 30 GeV 


5.7- 10~ 3 


1.5 • 10~ 2 


1.0 • 10~ 2 


2.3- 10~ 2 


8.6 • 10" 2 


cos(A(j) l T -jet) > -0.9 


4.7 • 10~ 3 


1.4- 10~ 2 


9.0 • 10~ 3 


2.2 • 10~ 2 


8.5 • 10" 2 


Re r-jet < 2.8 


4.5 • 10~ 3 


1.4- 10~ 2 


8.8 • 10~ 3 


2.1 • 10~ 2 


8.5 • 10~ 2 



Table 7: The cumulative acceptance in the mass window rnn ± 20 GeV, for different 
approaches of modeling production process. 
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7 Conclusions 



In the paper we have discussed one of the key signatures at LHC: (£ £p™ lss ) and (£ r— jet pjp iss ). 
They are important for searches of the Higgs boson if it decays into r-lepton pair. Nu- 
merical results were collected in sections 5 and 6. 

We have concentrated on the discussion of the reconstruction efficiency and mass res- 
olution of the LHC experiments. We have varied the production processes, also, we have 
used different methods for the simulation of the same channels, effectively re-shuffling 
parts of the hard process into parton shower. In this way we have obtained different 
event topologies. Even though we have applied rather simplified analysis, for example we 
have neglected additional cut-offs related to the background suppression and we have used 
simplified simulation for the detector response only, our analysis give an insight into the 
complexity of the problem. We have shown that the choice of production hard process 
affects the cumulative acceptance in the mass window for the signal reconstruction by a 
factor of few (even up to a factor of 10 for VBF fusion channel). The main effect comes 
from interplay of average transverse momentum of the Higgs boson, simulated differently 
depending on the way how the hard process is chosen and ambiguities in measurement of 
p-mtss _ ^his have an impact on the efficiency for finding physical solution of the neutrino 
system, necessary for the m TT reconstruction. Transverse momentum determines quality 
of the resolution, in particular in the regions of tails. 

Strong sensitivity on the production topology indicates that for the more complicated 
production processes, such as bb — > H Yukawa induced mechanism, the possible large 
theoretical systematic error may still need to be assumed for the predictions, even if 
small theoretical uncertainties obtained for the overall normalisation can be nowdays well 
controlled with NNLO calculations. For inclusive quantities effects related to the QCD 
regularisation and renormalisation scale are well understood. However, the fully exclusive 
Monte Carlo implementation is mandatory for the experimentally useful analysis. To 
convince the reader, we have explicitely studied effects of principal selection cut-offs which 
experiments will be most probably forced to use in data analysis. 

Let us stress, that we have not even touched the subject of the impact of those selection 
criteria which must be applied to optimise rejection against expected background (b-jet 
identification or b-jet veto). It will only add complexity to the discussion. Asking for b-jet 
identification or b-jet veto will enhance contribution from some NLO or NNLO terms, 
with respect to the total inclusive cross-section, in a manner which is difficult to predict 
analytically, as predictions need to be convoluted with detector effects (jets reconstruction 
and jets identification). Thus even stronger argument for the availabity of the existing 
theoretical calculations in the form of exclusive Monte Carlo, suitable for more educated 
studies than the one presented here. Without such a tool it will be difficult to decide 
in unambiguous way, how much of the mentioned above, factor of 3 to 10, ambiguity 
on cumulative acceptance, indeed contributes to theoretical uncertainty for the realistic 
signature of the Higgs boson. Note, that ambiguity depends drastically on how complete 
is the choice of selection cuts. If we had taken only first step (first lines of tables 2 and 
5) of the basic selection cuts, then we could conclude that acceptance is independent, 
both from different approaches to modelling of production process and from production 
process as well. Thus that it can be represented by an universal factor, up to about ±10% 
precision level. The conclusion would be, as we could see, overhelmingly false. On the 
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other hand, turning the question around, one can ask, if experimental selections can be 
modified to diminish dependence on the topological feautures of the events (namely quality 
of reconstruction and absolute value of the pjp iss ). 

Let us note also, that for the Higgs boson at 120 GeV dominant background will 
come from the irreducible Drell-Yan Z/j* — > rr production. The signal will appear 
above steeply failing edge of the resonant peak at the nominal Z-boson mass, and will be 
dominated by the misreconstructed on-shell Z-boson events. Calibration of the background 
from the shape outside signal mass window will be extremely difficult. 

Finally, let us recall again, that for the MSSM Higss scenarios in the parameter space 
corresponding to the A/H/h masses in range of 100-200 GeV, one will have to combine 
distinct production modes (gg — > H and bb — > H), and different final states of h — ► 
t + t~ with subsequent decays either into lepton-lepton or lepton-hadron channels. The 
signatures of the Higgs bosons h/H/A may overlap (even if there is no mass degeneracy 
between h/H/A). The goal of the experiment will be not only to establish evidence for the 
signal, but also to measure properties of the model (eg. couplings) from the observed versus 
predicted signal rates. The key challenge for evaluating remaining theoretical systematic 
error will require not only the proper normalisation of total cross section but also modeling 
of the production mechanism including differential distributions, essential as demonstrated 
here, for experimental signal reconstruction. 
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